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estimates of shear viscosity corrections at freeze-out, which we find important only for the elliptic 
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Figure 1: Distribution of charged particles in pseudorapidity at different centralities, calculated from the 
3+ ID hydrodynamic model compared to experimental data of BRAHMS and PHOBOS Collaborations. 

The expansion of the dense and hot fireball created in relativistic heavy-ion collisions can be 
described using relativistic hydrodynamics We perform 3+ ID hydrodynamic simulations of 
Au-Au collisions at the highest RHIC energy y/s = 200GeV Compared to existing calculations 
[^], we use a very short initial time 0.25fm/c and a realistic equation of state of dense matter 
The absence of a soft point in the equation of state leads to a rapid expansion of the system, resulting 
in a strong build up of the transverse flow. The initial energy density profile in the transverse plane 
is taken from the Glauber Model, and the initial distribution in space-time rapidity is adjusted to 
reproduce the measured charged particle distributions (Fig. [T]). 

Hydrodynamic equations of a perfect fluid d^T^^ = are solved for each impact parameter. 
The shape of the constant temperature {Tp = 150MeV) freeze-out hypersurfaces, as well as the 
final collective velocities of the fluid are exported to a statistical emission and resonance decay code 
THERMINATOR [|]. The THERMINATOR code generates complete events including statistically 
emitted particles from the fireball. Calculated transverse momentum spectra of produced particles 
follow very well the experimental results [||]. It is true for particle emitted at central rapidities for 
different collision centralities, up to 50%; also the agreement with the data at forward/backward 
rapidities is striking (solid lines in the left panel of Fig. ^ and proves that the hydrodynamic 
expansion model of the fireball combined with statistical emission applies in a broad range of 
rapidities. Microscopically it means that approximate local equilibrium is reached and maintained 
during a sizable time-span of the collective expansion. HBT correlation radii [§] calculated from 
the model are within 10% from the measured values (right panel of Fig. |^. A short initial time 
of the evolution requires a narrow initial distribution in space-time rapidity, this results in a strong 
suppression of elliptic flow when approaching the fragmentation regions (Fig 

We estimate the effects of viscous corrections on particle emission. Using the velocity flow 



2 



Flow, spectra and HBT radii in heavy-ion collisions 



Piotr Bozek 



10= 



> 

<u 
(5 

>; 

% 10-^ 

I- 
Q. 

CM 



■o 10' 



10-^ 



10 



10 



10 



10 



■13 



16 




4jr' 4ji 



\ I I I I I I L 

0.5 1 1.5 2 2.5 



[GeV/c] 



3 

o 



6 
4 
6 



0) 



c 
o 




1.5 



0) 
'<0 



3 

o 



0.5 



E 10 - 



1 - 



J 




c=0-5% 


1 


1 


* 

, 1 , 1 


■ 


■ 


, 1 , 1 






3+1Dti/s=0.08 

--■ 3+1Dti=0 
3+1Dti/s=0.16 


1,1,1,1 






— ^ ^. 


■ 

1 


a 

1 


* ■ 

, 1 , 1 



0.2 0.3 0.4 0.5 

[GeV/c] 



Figure 2: (Left panel) Transverse momentum spectra of at different rapidities for three different shear 
viscosity coefficients at freeze-out (r\ / s — Q,-^, j^, solid, dashed, and dotted lines) compared to BRAHMS 
Collaboration data. (Right panel) HBT correlation radii for different strengths of shear viscosity corrections 
at freeze-out, STAR Collaboration data. 



obtained in a 3 + ID perfect fluid dynamics we calculate nonequilibrium corrections to the Cooper- 
Frye formula from shear viscosity [^. Modified momentum distributions are implemented in the 
statistical emission code. Results range from, no appreciable modifications of the transverse mo- 
mentum spectra for p± < 1.5GeV and central rapidities, up to a 30% increase at rapidity y = 3.5 
and p± = 3GeV for rj/s = 0.16. The HBT radii are not sensitive at all to the shear corrections at 
freeze-out if the the flow remains unchanged (Right panel in Fig. ^. 

A significant reduction of the elliptic flow is induced by stress corrections (Fig. |3|). The 
reduction is 20% at central rapidity and goes up 60% at pseudorapidity 4 for T]/^' = 0.16. This 
observation agrees with the results of Ref. where dissipation from hadronic rescattering was 
found to be important at large rapidities. This effect modifies strongly the dependence of the elliptic 
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Figure 3: (Left panel) Elliptic flow coefficient for charged particles as function of pseudorapidity for three 
different viscosity coefficients at freeze-out, compared to PHOBOS Collaboration data. (Right panel) Ratio 
of the stress velocity gradient to the overall expansion rate as function of time at two space-time rapidities. 
The arrows indicate the time of the freeze-out. 

flow of charged particles on pseudorapidity. From the ratio = — presented in 

the right panel of Fig. ^ we expect large shear viscosity corrections at large rapidities, whereas 
at central rapidities correction from bulk viscosity (proportional to dau"') could also be important. 
The approximate agreement of perfect fluid calculations with the data on elliptic flow is accidental. 
Moreover shear viscosity modifies also the longitudinal acceleration of matter changing the 
distributions in Figs. ^ and ^ 

Summarizing, we find good agreement of the results of 3 + ID perfect fluid hydrodynamics 
with the measured transverse momentum spectra, HBT radii and elliptic flow. Estimates of shear 
viscosity effects show that the last observable is not robust and a reliable estimate thereof requires 
the use of a fluid expansion model including shear and bulk viscosity. 
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